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Abstract 

The spatial behavior of an animal population of dispersive and competitive species 
is studied. It may be realized as a result of various kinds of biological effects, as hetero- 
geneity of environmental conditions, mutual attractive or repulsive interactions between 
individuals, competition between individuals for resources and localization of the off- 
spring during a reproduction event. A stochastic model describing the behavior of each 
individual of such a population is introduced and a large population limit is studied. As 
consequence, the global existence of a nonnegative weak solution to a multidimensional 
parabolic strongly coupled model for competing species is proved. The main feature of 
the model is the nonlocal nonlinearity appearing in the diffusion and drift term. The 
diffusion matrix is non-symmetric and generally not positive definite and the cross- 
diffusions terms are allowed to be large if subspecies sizes are large. We prove existence 
and uniqueness of the finite measure- valued solution and give assumptions under which 
the solution takes values in a functional space. Then we make the competition kernel 
converge to a Dirac measure and obtain solution to a locally competitive version of the 
previous equation. The proofs are essentially based on the underlying stochastic flow 
related to the dispersive part of the dynamics. 



1 Introduction 

The spatial structure of a biological community is a fundamental subject in mathematical 
ecology and in particular the spatial distribution formed by dispersive motions of pop- 
ulations with intra- and inter- specific interactions. (See for instance Keller-Segel 
Gurney-Nisbet [9], |10j . Mimura-Murray |17| . Shigesada-Kawasaki-Teramoto [18J). In this 
paper, the spatial behavior of a population of competitive species is studied. The disper- 
sive motion of an individual in its environment is realized as the result of various kinds of 
biological effects, as heterogeneity of environmental conditions, mutual attractive or repul- 
sive interactions with the other individuals, competition for resources and localization of 
the offspring during a reproduction event. These different effects will be modeled by local 
interaction kernels depending on the type of the individual and acting either on its spatial 
parameters or on its ecological parameters. 
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The population is composed of M sub-populations (species) characterized by different phe- 
notypes. Each species has its own spatial and ecological dynamics depending on the spatial 
and genetic characteristics of the whole population. We assume that each individual (of a 
given type) motion is driven by a diffusion process on M. d whose coefficients depend on the 
spatial repartition of the different species around. Moreover, the individuals may reproduce 
and die, either from their natural death or because of the competition pressure for sharing 
resources. Each species has its own growth rate. The competition pressure of an individual 
of type j on an individual of type i depends both on the location of these individuals and 
on their type. It is not assumed to be symmetric on i,j. If i = j, it is called infra-specific 
competition and inter-specific in case i ^ j. 

In the paper, we describe the stochastic dynamics of such a population by an individual- 
based model. Each individual is characterized by its type and its spatial location. Because of 
the births and deaths of individuals, the population doesn't live in a vector space of positions 
and we model its dynamics as a Markov process with values in the M-dimensional vector 
space of M d -point measures. We introduce the charge capacity parameter K describing 
the order of the population size. To be consistent, the individuals are weighted by 
The existence of the population process is obtained by standard arguments. Then large 
population asymptotics is studied. We show that when K tends to infinity, the population 
process converges to a weak solution of the nonlocal (in trait and space) nonlinear parabolic 
cross-diffusion-reaction system defined by: for all i € {1, • • • , M}, 

d 
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Here, u l (t, .) denotes in general a finite measure on R" for any t > 0, and {G^ ,H ij , C' J ')i<i,i<M 
re 3M 2 -nonnegative and smooth L 1 -functions defined from IR d to R+ that model the spatial 
interactions between individuals of type i and j. They are not assumed to be symmetric. 
By this convergence result, we prove a theorem of existence of a weak solution for Equation 
dU) . Next we prove the uniqueness of such a solution and we give two sets of assumptions 
for which the measure solution has a density with respect to the Lebesgue measure, thus 
ensuring the existence of a function solution to ([I]). The tools we use are probabilistic ones 
and rely on the use of the stochastic flow related to the dispersive part of the equation and 
its inverse (as function of the initial position). 

In the model leading to Equation ([1]), the competition between two individuals is described 
as a function of the distance between the individuals. This biological assumption is clear: 
the closer the animals are, the stronger is the fight to share resources. An extreme situation 
is the local case for which individuals only compete if they stay at the same place. Math- 
ematically speaking, it means that for i,j £ {1, • • • , M} and for x 6 K rf , the competition 
kernel has the form 

C i i(x) = c ij C £ (x), 

where are positive constant numbers and the measures C e (x — y)dy weakly converge to 
the Dirac measure at x when the range of interaction e tends to 0. For instance, C e may be 
the centered Gaussian density with variance e. We study the convergence of the solution 
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u £ of ([T]) when e tends to zero. We show that u £ converges to the solution u of the spatially 
nonlocal nonlinear cross-diffusion equation: for all i S {1, • • • , M}, 

d 

dtu' = ^4(aij(.,G ll * tt 
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Here again, the tools are based on the stochastic flow of the underlying diffusion process. 

To our knowledge, the nonlocal nonlinear Equations (pQ) and (|2]) have never been studied, 
despite the fact that they naturally arise from the biological motivation. In Champagnat- 
Meleard [3], Arnold-Desvillettes-Prevost [1], Bouin et al. [2], models with space displace- 
ments and a continuous trait are studied, but the diffusion coefficients in these models only 
depend on the type of each individual, and not on the spatial distribution of the other 
animals alive. The dependance of the diffusion coefficient on the individual density is in- 
deed the main new difficulty we deal with in this paper and it requires the introduction of 
different techniques from those developed in the aforementioned works. 
Cross-diffusions with local spatial interaction and local competition have excited the scien- 
tific community, see for example Mimura- Murray [17] , Mimura-Kawasaki [16] , Lou-Martinez- 
Li [15], Chen-Jungel [5], [6] and more recently Lepoutre-Pierre-Rolland p3], Desvillettes- 
Lepoutre-Moussa [7J and references therein. The prototypical equation in this situation is 
given by 

dtu 1 =d\A ((ai + bi 2 u 2 )u l ) + (ri - cnu 1 - ci 2 n 2 )u 1 , 
d t u 2 =d 2 A ((o 2 + b 2 iu 2 )u 1 ) + (r 2 - c 2 \u l - c 22 u 2 )u 2 , 

with boundary conditions on a given bounded smooth domain of R . Although global 
existence results were obtained in [5], [6] for such equation, their techniques do not seem 
to be useful in more general situations. In that direction, an interesting but highly difficult 
open challenge is to establish the convergence of the system ([2]) to systems of the type of 
([3|) when the kernels G and H tend to Dirac measures. 

2 The Individual-based Model 

2.1 Assumptions 

For i = 1, • • • ,M, let us denote by S + (R d ) the space of symmetric nonnegative diffusion 
matrices and define the measurable functions 

a* : R d x R M ^ S + (R d ); 
V : R d x R M ^ R d ; 
ri : R d ^ E + . 

We denote by a 1 the d x d— matrix such that a 1 = a 1 (<r 1 )* . We will assume that a % : 
R d x R M ^ S + (R d ) and b { : R d x R M -> R d are measurable and satisfy 
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i) There is a positive constant L such that for any x, x' G M. d and any Vj,v'j, G R+, (j G 
|cr J (x,wi, • • • ,v M ) ~ cr l (x',v' 1 , ■ ■ ■ ,v' M )\ + |6 l (x,wi, • • • ,v M ) - b l (x',v' 1 , ■ ■ ■ ,v' M )\ 

M 

<L(\x-x'\+Y,\^-^\- 

3=1 



ii) There exists a positive constant Cm such that for any x£l' ! and any v\ , • • • , «m G M+ 
|cr(x,fi, • • • ,v M )\ < C M \x\. 

hi) The functions (G 1 ^ , H 1J ,C lJ )i<ij<M defined from M. d to M. + are assumed to be non- 
negative, continuous and bounded. 

iv) The nonnegative functions r\ are assumed to be bounded . An upper bound is denoted 
by fj. 



2.2 The M-type and diffusive stochastic population dynamics 

Let us now describe the dynamics of the population we are interested in. We take into 
account the births and deaths of all individuals and their motion during their life. 
The population dynamics will be modeled by a point measure-valued Markov processes. 
Let us fix the charge capacity and define 



as the space of weighted point measures on M. d . The stochastic population process (i't < )t>o 
will take values in {Mk) M ■ The ith coordinate of this process describes the spatial config- 
uration of the subpopulation of type i. Thus 
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For any i G {1, ■ ■ ■ , M}, N£ G N stands for the number of living individuals of type i at 
time t and x]' 1 , ,X^ 1 ' 1 describe their position (in R d ). 

The dynamics of the population can be roughly summarized as follows: 

The initial population is characterized by the measures (^o)i<«<m G {M.k) M at time t = 0. 
Any individual located at x G M. d with type i at time t has two independent exponential 
clocks: a reproduction clock with parameter ri(x) and a mortality clock with parameter 
*v 3 t ' (x). If the "clonal reproduction" clock of an individual rings, then it produces 
at the same location an individual with the same type as itself. If the "mortality" clock 
of an individual rings, then this individual disappears. The death rate of a ith individual 
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depends on the positions of the other individuals through the kernel C %3 which describes 
how species j acts on i in the competition for resources. 

During its life, an individual will move as a diffusion process whose coefficients depend on 
all individual positions. The motion of an individual with type i is a diffusion process with 
diffusion matrix a* (., G il * v l > K , ■■■ , G iM * v M > K ) and drift vector 6* (., /P 1 * v l > K , • • • , H iM * 
v M ' K ). The coefficients take into account the effects due to the nonhomogeneous spatial 
densities of the different species. Indeed, a species can be attracted or repulsed by the 
other ones and the concentration of species may increase or decrease the fluctuations in the 
dynamics. 

The measure-valued process (yt)t>o is a Markov process which can be rigorously expressed 
as solution of a stochastic differential equation driven by d-dimensional Brownian motions 
(B n ' l )i<i<M,neN* and by a Poisson point measure Q(dt,dn,dO) on W + x N* x M + with 
intensity dt Yl n m* ^ n ^> an independent and independent of the initial condition uq. The 
measure Q stochastically dominates the jump process which describes the births and deaths 
in the population. The Brownian motions drive the spatial behavior of the individuals alive. 
By using Ito's formula and for C 2 (M d )-functions fi, i = 1, • ■ ■ , M, we get 

1 K(vl,l) 

w K Ji) = ^ E W) 

n=l 

= (4, fi) + Jrj EE W * W)- • • • , G iM * W))) 
J ° n=l k,l 



d x J t (X^)dBr(s) 

+^ jf ("i. ( E ^ * ••• , G m * <)t,.,,.r,ixr) 

+ E b%(, H 11 * vl> K , ■■■,H iM * v?> K ) d Xk fi(X^)\ )ds 

+ / ,t]xNx[0,i] ( lo ^^ r 'W-) ~ 1 ^(^-)<^<^(^r-)+E J A i 1 ^^^(x^)) Q(ds,dn,d6). 

(4) 

The law of (r^)t is characterized by its infinitesimal generator L which captures the dy- 
namics described above. L is the sum of a birth and death (ecological) part L e and a 
diffusion part L^. The generator L e is defined for bounded and measurable functions <fi 
from (Mk) M into R. Let us denote by e« the ith unit vector of R M . For i/ = (i^)i<j<M = 
En=i <W)i<i<M, we define 

M N* ^ 

L e <p{ v ) = Y^r i {x n > i ) E + k$x^ -ei) - 0(f)) 

i=l n=l 



M N 1 M 

+EE (E^' * ^"(* B,< ))(^ - - (5) 
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To define the diffusion part of the generator we need to introduce a standard class of 
cylindrical functions generating the set of bounded and measurable functions from (Mk) M 
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into R. Let us consider F € Cf (M A/ ) and for i £ {1, • • • ,M}, we introduce f { £ C 2 
and define 

F / (i/)=F(( I / 1 , /!),-.. > M ,/m». (6) 

The diffusive part of the generator can easily be deduced from Ito's formula. Its form 
is similar to the one obtained in the whole space for branching diffusing processes and is 
given by 

M / d 



L d F f (u) = ^(j^(v\{j kl {-,C?Uv\...,G iM *v M )dl kX ^ 
i=i \k,l=i ^ ^ 

+ J2bU;H il *u\... ,H m *v M )d x jX)Fl{{v\h),--- ,{v M J M )) 
k=i J ' 

+ E ali(;G a *v\--- , G iM * v M ) d Xl fi d Xk fj\ F^v 1 , /i), • • • > M /a/»Y 



J 
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(7) 



Hence, 



LFf{v) = L e F f {v)+L d F f {v). (8) 



3 Large population approximation and non local Lotka-Volterra 
cross diffusion system 

3.1 Existence and uniqueness of weak measure solutions 

We now state a large population approximation for the previous M species model by making 
the charge capacity K tend to infinity. This result implies in particular the existence of weak 
solutions to a non local cross-diffusion system of nonlinear partial differential equations. 
We denote by M the space of finite measures in M. d endowed the the weak topology. 

Theorem 3.1 Assume that for somep > 2, sup^ ~E((i>q K , l) p ) < +oo for any i = 1, • • • , M. 
Assume also that the coefficients satisfy the assumptions of Section \2. 1\ and moreover that 
the sequence of finite measures (y^ , ••• ,v ' ) converges in law as K goes to infinity to 
the deterministic measures (£q,-"" >£o f ) • Then, when K tends to infinity, the sequence 
[y 1,K i ■ ■ ■ ,v M,K ) converges in law in O([0, T],A4 M ) to the unique deterministic continu- 
ous finite measure-valued function £ = • • • , £ M ) weak solution of the following cross- 
diffusion system: for all i = 1, • • • , M and all C^' 2 -function f % = fl{x), 

<& ft) = fo) + [ J { E a «(-. * & • • • > ^ * ^) d l*Js 

M . 

+ £ bU; H* 1 H iM * ef) d Xk f s + (n - E * &tt + dsf s (x) £(dx)ds. 

k j=l J 

(9) 



6 



Remark 3.2 Observe that in contrast to the models of cross-diffusion introduced in Jiff] / or 
fT(j\/ , Equation © allows for long range interaction in the coefficients of spatial diffusion. 
For example, taking G tJ = H 13 = 1, the spatial behavior of individuals of type i depend 
on the total mass of the subspecies j. It also covers some cases where the diffusion matrix 
might vanish , e.g. a l (x, v\, . . . , vm) = Id^i(52jLi v j) with ^ : [0, oo] — > R+ a smooth 
bounded function vanishing at and Id the identity matrix. 

Remark 3.3 Notice that a solution to satisfies sup ig [ 0T ] 1 1 CI 1 1 TV < eriT ||£ollTV f or 
i = 1, M as is readily seen by taking f % = 1 and using the non negativity of the functions 
C lJ and Gronwall's lemma. 



The proof is decomposed in several steps: propagation of moments, uniform tightness of 
the laws of [iA> , • • • , v M,K ), identification of the limit and uniqueness of the solution of 
([9]). The three first steps are standard and can be adapted from the arguments in Fournier- 
Meleard [8] and Champagnat-Meleard [3]. However, the proof of uniqueness requires new 
arguments to deal with the interaction in the diffusion terms. That result as well as the 
convergence of the previous model to a local competition model will rely on the use of 
stochastic flows which we now recall. 

Given a solution , . . . , £t)te[o,T\ to ([SJ, we set for % = 1, . . . , M, 

a(i,t,x) = a\x,G a *CHx),--- ,G ai *^(x)) 

b(i,t,x) = VfaH**^),--- ,H iM *£* s (x)). (10) 

We then consider for each x G M. d and s G [0, T] the unique (strong) solution 

Xl t (x) = (Xlj(x),...,Xl> d t (x)) te[s,T] 

of the stochastic differential equation 

Xi >t (x) = x + ! a{i,r,Xi r (x))dB l r + ! b{i,r, Xi >r (x))dr, t G [s,T] (11) 

J s J s 

where B l = {B l ' q )^ 1 is a standard cf-dimensional Brownian motion in a given probabil- 
ity space. The three parameter process (s,t,x) \-t X l st (x) is called the stochastic flow 
associated with the coefficients o~(i,t,x) and b(i,t,x). 

Uniqueness of weak measure solutions of ([9]) is established in the next proposition. For 
notational simplicity, we will deal only with the case M = 2. All arguments easily extend 
to the general case. In all the sequel C > will denote a constant that may change form 
line to line. 



Proposition 3.4 Let £ = (C 1 ,^ 2 ) and £ = (C 1 ^ 2 ) be two solutions of the system ([9]) in 
[0,T] with M = 2. Then = for all t G [0,T]. 

To prove Proposition 13.41 we shall use a norm which is weaker than to total variation one, 
yet better adapted to perturbations in the diffusion coefficients than the latter. We denote 
by CB(R. d ) the space of Lipschitz continuous and bounded functions on M d , and by || • \\cb 
or simply || • || the corresponding norm, 

\W\\CB ^sup '^-y + su P K*)|. 

x^y \ x y\ x 
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We endow A4(M. d ) with the dual norm with respect to CB(M. d ), 

\\r)\\cB* ■= sup (ri,cp), J? £ M^) 

Given two solutions £t)te[0,T\ and (£h£i)te[a,T\ of &, let ^,t0) and -X^.tO 15 ) be the 
processes defined in (fTTj) and define analogously coefficients cr(i, i, x) and b(i, i, x) in terms 
°f (? 1 )^ 2 )) an d the stochastic flows X l st {x) given for i = 1,2 by the solution to the SDEs 

X l st {x)=x+ f a(i,r,X i s , r (x))dB t r + [ b(i,r, X i 8<r (x))dr, 

J s J s 

respectively driven by the same Brownian motions B 1 and B 2 as Xj t (x) and X 2 t (x). 
We also need some facts about the non homogenous transition semigroups of X\ t (x) and 
X\ t (x), which we respectively denote by P l st {x,dy) and Pg t (x,dy). 

Lemma 3.5 For allT>0, ipe CB(R d ), < s < t < T , i = 1,2 

a) \\Pi,t<P\\cB < C(e,e,T)\\cp\\l B and ||P^|| 2 < C(£V 2 ,T)IMI£b . 

b) For all x G R d , 

\pIMx) - PlM x )\ 2 < c'(t\t\t\~t\T)y\\ 2 CB f - £l\\l B * + \\g - ~g\\ 2 CB *dr 

J s 

The constants depend on^ 1 (resp. t^ 1 ) only through sup te r 0T i IICIIItv (resp. sup tg [ ,T] II£!IItvJ- 

Proof: a) Is is enough to control the Lipschitz constants of the function P£ t <p or P % st ty- 
We have, by Burkholder-Davis-Gundy inequality 



e\KM - KM\ ^°\ x - y\ + c / r ' X IA X )) - r > x l,M)\ dr 

J s 

+ C I E\b(i, r, XUx)) - b(i, r, X*Jy))\ 3 dr 



<C\x-y\ z + C / E\Xl >r (x) - Xl jr (y)\ dr 

J s 

for all s < t and constants depending on the Lipschitz constants of the coefficients a 1 , b 1 
and of the functions G 3 and H lj . By Gronwall's lemma, t (x) — X\ t (y)\ 2 < C\x — y\ 2 
which easily yields 

\PlMx)-PlMy)\ 2 <cM 2 CB \x- y \ 2 

as required. 

b) For notational simplicity we consider first the case 6 = 0. Using similar types of inequal- 
ities as before, we have for all s < t < T, 

E\Xt >t (x) -Xl t (x)\ 2 < C'(j\\Xl r (x) -Xl, r {x)\ 2 + E\G il *6PC») - & 1 *UKA X ))\ 

+ E\G* 2 * UXIA X )) ~ Qt2 * i2(Xl r (x))\ 2 dr^j 
<C"(j\\X\ r {x) - Xl r {x)\ 2 + E J G a (Xi r ^)-y)(e s (dy)-iHdy)) 
+ JE / G a (XiJx)-y)(e s {dy)-l\dy)) <lr 
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Since the functions y t- > G lj {X l s (x) — y) are uniformly Lipschitz continuous, we deduce with 
Gronwall's lemma that 



E(\Xl tt (x) - Xl t (x)\') <C / HC 1 - e\\ 2 CB* + IIC 2 - eXcB*dr 

J s 

which allows us to easily conclude. The case b 7^ is similar with additional terms involving 
the kernels H l K [] 



Proof of Proposition 13. 4t Take (p £ CB(R ) with \\(p\\ < 1. By the Feynmann-Kac 
formula, the function f^(s,x) = ~E(<p(Xl t (x))) = Pg t <p(x) is the unique solution of the 
linear parabolic problem 

d s fV(s, x) + aU; G 11 * e t , G 12 * g)d XkX J {t \s,x) + bl(-,H 11 *e t ,H l2 *g)d Xk / t \s,x) = 

with final condition at time s = t equal to f(x). Replacing /w in the first equation in ([U]), 
we see that £ x satisfies 

<6\ <P) = {tlKtV) + f I Hx) - C u * e s (x) - C 12 * e s (x))PlMx)£(dx)ds, (12) 
and, similarly, 



(#,¥>> = P » + J I (ri(x) - C n * H(x) - C 12 * i 2 (x))PlM^H(dx)ds. 
Consequently, 



(Pl t -Pl t )ip(x)e s (dx 
2 



+ 



+ 



c n *(e s -ii)(x)PiM^e s (dx 

C 11 * H(x)PlMx)(e s (dx) - H(dx)) 

C 12 * i 2 s (x)(Pl Mx) - PiMmlidx 



PlMx)(e s (dx)-g(dx)) 

c 11 * il(x)(PlMx) - PlMx))e s (dx) 2 

c i2 *(e s -£ 2 s )(x)p}Mx)e s (dx) 



+ 



C 12 * i 2 (x)P}Mx)(e s (dx) - H(dx)) 



<C sup \(Pl t -Pl t My)\ 2 + C / snp\(Pl t -Pl t My)\ 2 + (e s -il,Pl^) 2 
y JO I 9 

+ sup|C 11 *(£ s 1 -£)(y)| 2 +' ' 

V 



+ sup|C"*(^-^)(l/)r + 
y 



C 11 * H(x)PlMx)(£(dx) - H(dx)) 
C 12 * £ 2 (x)PlMx)(e s (dx) - H(dx)) 



-1 2 



ds, 



for constants depending on sup te r 0)T ] ||C?IItv su Pte[o,T] ll£tlllv an< ^ T- ^ ne functions x \-± 
C l i(x—y) are Lipschitz continuous uniformly in y and, by Lemma f3.5l a). C 1 ^ *^i{x)P} t ^p{x) 
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are Lipschitz continuous bounded functions, uniformly in s,t G [0, T]. Together with 
Lemma 13.51 b), this entails 




and we can analogously obtain a similar bound for (^ 2 — ^ 2 , <p) 2 ■ Taking sup 1 1 y> 1 1 < 1 > summing 
the two obtained inequalities and using Gronwall's lemma we conclude that 

\\Ct-Ct\\a^ = 

for all t £ [0, T] and i = 1, 2 and thus uniqueness for System ([9]). [] 
3.2 Regularity of the flow and function solutions 

We next show under two types of suitable assumptions on the coefficients and the initial 
condition that for i = 1, ...,M the solution Q{dx) has a density Q(x) with respect to 
Lebesgue measure. 

Lemma 3.6 Let . . . , £i^)te[o,T] ^ e measure solution of Q. Then, for all t G [0, T], 
where m\ is defined by 

(ml if) := E (7 ^ )t (*))^)) (13) 
/or any bounded function ip. 

Proof: We write the proof for M = 2 and * = 1, and omit for notational simplicity the 
superscript 1 in the flow X\ t {_x). Taking in the first equation in ([9]) the function 

f®(s,x) :=E(<p(X a , t (x))explJ ri(X a , r (x)) - C n * £ {X Sir (x)) - C 12 * ^{X s>r (x))drV\ , 

which (by the Feynman-Kac formula) is the unique solution of the parabolic problem 

=dsf®(s, x) + G il * G a * e t )d, kX J (t) (s,x) + bi(-,H 11 * e t , H l2 *Ht)d x J (t) (s,x) 

+ ( ri (x)-c ii *es-c i2 *es)f {t) (s,x) 

with final condition at time s = t equal to p(x), we get that 

{Zl<p) = E (J d <p(X , t (x)) exp { jf niXo^x)) - C u * &(X , r (x)) - C 12 * g(X 0>r (x))dr\ $ 

(14) 

for each continuous bounded function p > 0, so that ip) < e riT (m\, ip) for all bounded 
continuous cp. The measure + e riT m\ being regular, for each Borel set A and e > we 
can find a closed set B C A s.t. — e riT m\,A\B) < e. Since the sequence of bounded 
continuous functions fk{x) = (1 — kd(x,B)) V pointwise converges to 1b as k — > 00, 
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(£t — e riT m\,\B ~ fk) by dominated convergence w.r.t. the positive finite measure 
|£j — e riT m\\ . It follows that for k sufficiently large 

(& 1 - e^ T ml,A) < + m],A\B) + (& 1 - e^mj, 1 B - f k ) + (& 1 - e^ T m\j k ) < 2s, 

that is, A) < e flT (m t 1 , A). □ 

We immediately deduce 

Corollary 3.7 For any initial measure condition (£q, • • • , £oO, * n ^ e uniform elliptic case: 
3\i > suc/i i/toi y*a l (x,v)y > \{\y\ 2 Vx,y G K^f G i/ie measure Q(dx) has a density 
£l(x) with respect to Lebesgue measure for all t G (0, T]. 

As pointed out in Remark 13.21 some natural biological examples are not covered by this 
ellipticity assumption. We provide a finer result covering some non elliptic cases under 
additional regularity assumptions. 

Let us denote by C^' a the space of twice differentiable functions with an a-Holderian second 
derivative (and a similar definition for C^' a ). 

Proposition 3.8 Assume for each i = 1, . . . , M that o~ l (x, v\, . . . , vm) and b l (x, v%, . . . , vm) 
are respectively globally C^' a (M. d x [0, oo) M ) and globally C^' a (lR d x [0,oo) M ) for some 
a G (0,1). If for some type i the measure £q has a density, then Q has a density for 
all t G [0,T]. 

We will need some classical properties of the stochastic flow stated in Kunita [12] that we 
summarize in the next Lemma. 

Lemma 3.9 Under the assumptions of Proposition [J7E[ the process (s,t,x) t-t X l st {x) has 
a continuous version such that, a.s. for each s < t the function (s,t,x) t-t X l st {x) is a 
global diffeomorphism of class C 1 '^ for all ft G (0, a). 

Moreover, for each (t, y) G [0, T) x R d the inverse mappings rf s t (y) :— (X^ t ) 1 (y), < s < 
t <T, satisfy the stochastic differential equation 

Vlt(y) = y- [ v(i,r,4,t(y))dB l r - I b(i,r,4,t(y))dr (15) 

J s J s 

where b(i,r,y) = b(i,r,y) — Ylt i=i °~ik{h r > y)d yi o-,k{h T, y) and dB l refers to the backward 
ltd integral with respect to the Brownian motion B % (see Kunita [12] p. 194)- 

Finally, for each (t,y) G [0,T]xIR d the (invertible) Jacobian matrix V y rf s t (y) = (dy^J^yfj 
of rf st (y) satisfies the system of backward linear stochastic differential equations : 
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Proof: The assumptions on the coefficients a 1 and b l and on the kernels G y and H lJ , 
together with the fact that sup tg j ,T] I Id 1 1 TV < oo for i = 1,....,M imply that for each i the 
functions x i— > o~(i,t,x) and x i— > b(i,t,x) are respectively C 2 ' a (M. d ) and C 1,a (M. d ) for some 
a £ (0,1), uniformly in t E [0, T]. Consequently, we can apply Theorems 3.1 and 6.1 in 
Kunita [T2]. 

□ 

Proof of Proposition 13.8b Let ip > be a bounded measurable function in W d . By the 
previous lemma, we can do the change of variable Xq^(x) = y in the integral inside the 
expectation in (fl3|) to get: 



oo > (m t , (p) = E 



/ ( /'(y)Co( 7 ?o,t(y))|rfetV77o,t(?/)|dy 



Fubini's theorem (everything being positive) yields that mt has the (integrable) density 
y i— y E [Coivo t(y))\detVr]o t {y)\] with respect to Lebesgue measure. 

□ 



4 Convergence to local competition 

Our aim in this section is to describe some situations where the interaction range of the 
competition is much smaller than the one of spatial diffusion. For example one may assume 
that animals interact for sharing resources as they are on the same place but diffuse depend- 
ing on the densities of the different species staying around them in a larger neighbourhood. 
To model such situation, we suppose now that C lJ = c l3 ^ e for c lJ > some fixed constant 
and 7 e a suitable smooth approximation of the Dirac mass as e — > 0. Our goal is to show 
that, under additional regularity assumptions, the (unique) solution £ = • • • ,£ A/ ) of 
Equation ([9]) given by Theorem 13. II for such competition coefficients converges as e — > 0, to 
a weak function solution of the system of Equations ([2|) . 

Theorem 4.1 Assume that for each i = 1,...,M the coefficients satisfy the assump- 
tions of Section [2A\ and that o l {x, v±, . . . , vm) and b l (x, v\, . . . , vm) are respectively globally 
Cfr' a (M. d x [0,oo) A/ ) and globally C^' a (M. d x [0, oo) M ) for some a G (0, 1). Assume moreover 
that the functions r% have bounded derivatives, and that the initial measures (£q,--- >£o ) 
have bounded densities with bounded derivatives. 

Furthermore, assume that C* J = c l ^ e for c 4J > some fixed constant and 7 £ = ^{x/e)e~ d 
for some regular function 7 > satisfying f Rd ~/(x)dx = 1 and j„ d \x\^{x)dx < 00. 
Then, for each T > the unique weak function- solution £ e to Equation ([9]) converges in 
the space C([0,T], Ai M ) (endowed with the uniform topology) at speed e with respect to the 
dual Lipschitz norm to a solution u = (u 1 , • • • , u M ) of the following non local cross- diffusion 
system with local competition: 

(4, fl) = <& /5> + jT / { E °«(-' G il *u\,--- , G iM * uf ) d XkXl fl 

M . 

+ £ b{(; W 1 * u\, ■ ■ ■ , W M * uf 1 ) d x Jl +{n-Y, <-' J "i)-i: + dsfl Ux) u\{x)dxds. 

k j=l > 

(17) 
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Moreover the function u is the unique function solution of (|17h such that sup te r 0i T] [|itt[|i + 
\\ut\\ct3 < oo. 

We will extend techniques previously used in the uniqueness result to a convergence argu- 
ment using the same dual norm and some additional estimates and technical results. 

Lemma 4.2 Under the assumptions of Theorem \4-l\ for each t G [0, T] the functions 
(£f i = 1, . . . , M have bounded first order derivatives. Moreover, there exists for each i a 
constant K \ > depending on the functions C u ,j = 1, • • • , M only through their L 1 norms 
c lJ (and not on e), such that 

max{ sup ||(£Dilco, sup \\V(^Y IU} < K h i = 1, . . . M. 
te[o,T] te[o,T] 

The previous result relies on some refinements of Lemma 13.91 

Lemma 4.3 Under the assumptions of Lemma \3. Sh for each % = 1, . . . , M and p > 2, there 
exist finite constants Kn{p) > and Ki2{p) > not depending on the nonnegative kernels 
C ij such that for all t G [0,T], 

supE sup \V y r]l t (y)\ p ) < K a (p) and sup E sup \detV y r]ij(y)\ p ) < K i2 (p). 

y£R d \s€[0,t] / yeR d \s£[0,t] J 

If moreover the assumptions of Lemma \4-%\ are satisfied, then for each s < t with s,t £ 
[0,T] the function y \— > detSJ y r] l s t (y) is a.s. differentiable and there exists Ki^(p) > not 
depending on the nonnegative kernels C 13 such that for all t G [0, T], 

sup E ( sup \V y [detVyrii/y)] \A < K i3 (p). 

y£R d \s€[0,t] J 

Proof of Lemma 14. 3t For fixed i G {1,...,M} and t G [0, T] we define coefficients 
P : [0,t] xR^R and A q : [0, t] x R d -)■ R, q = 1, . . . , d with components @ k and A\, 
fe = l,...,dby 

P(s, V) ■= ~b(h t - s, y), A q (s, y) := a. q (i, t - s, y), 



and the process (Z s (y);s G [0,t],y G R ) by Z s (y) = i]l_ st (y). Then, denoting by W = 
(W 1 , . . . ,W d ) the standard d— dimensional Brownian motion W s := B\_ s — Bf, it easily 
follows from Lemma 13.91 that Z s {y) satisfies the classic ltd stochastic differential equation 



s rs 



Z s (y)=y+ A(r,Z r (y))dW r + (3(r, Z r {y))dr (18) 
J o J o 

whereas the associated Jacobian matrix satisfies the linear system: 

d yi Z k s (y) = I d + I* d Xp Al(r,Z r (y))d yi ZP(y)dW? + J2 f d Xp /3 k (r, Z r (y))d Vl Zf(y)dr. 

(19) 
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Using the Burkholder-Davis-Gundy inequality, the boundedness of the derivatives of A q 
and (3 and GronwalFs lemma, we deduce that 



E sup \V y Z s (y)\P < K a (p) 
\se[o,t] J 



(20) 



for some constant Kn (p) which depends on bounds for those derivatives and on sup ig [ ,T] llCIIIrVj 
and does not depend on y £ R^. This yields the first asserted estimate. 
In order to get the estimates for the determinant and its gradient, we rewrite (|18p in 
Stratonovich form 



Z s {y) = y+ [ S A(r,Z r (y))odW r + f 'p°(r,Z r 
Jo Jo 



(y))dr 



where /3°(r, x) = /3(r, x) — | Yli q ^f( r J x )9 yi A q (r, x). By the proof of Lemma 4.3.1 in Kunita 
[13], de£V y Z s (y) satisfies the linear Stratonovich stochastic differential equation 



detV y Z s (y) = 1+1 detV y Z r (y) \^ 
Jo k=i 



£ d yk Al(r, Z r {y)) o dW? + d y J° k (r, Z r (y))dr 

q=l 



(21) 

Using Burkholder-Davis-Gundy inequality in the ltd form of the previous equation, we 
deduce using also Gronwall's lemma that 



E sup \detV y Z s (y)\ p < K i2 (p) 
\se[o,t] J 



(22) 



for some constant Ki2(p) depending on bounds on the (up to third order) derivatives of a 1 
and (up to second order derivatives) of b and on sup tg r 0T ] ||£!IItv; this yields the second 
required estimate. Remark 13.31 ensures that under the assumptions of Theorem 14. 1\ the 
constants Kn and K{2(p) don't depend on e. 

Finally, under the assumptions of Lemma 14.21 we deduce from equation (|2ip and Theorem 
3.3.3 of [13] (see also Exercise 3.1.5 therein) the a.s. differentiability of the mapping y \— )■ 
defVyT/l t (y), and that its derivative with respect to yi satisfies 



d yi [detV y Z s (y)}= / d yi [detV y Z r (y)]Y^ 



+ f S detV y Z r (y) £ d m ZT(y) 
Jo 



m,k=l 



k=l 
d 



d yk Al(r, Z r (y)) o dW? + d Vk ft(r, Z r (y)) 



9=1 



£ d 2 ymyk Al(r, Z r {y)) o dW< + d ym yj° k (r, Z r {y)) 

9=1 



dr 



dr. 



Writing this equation in ltd form, we now deduce with the Burkholder-Davis-Gundy in- 
equality that 4>(s) := E (sup rg p,sl [detV y Z r {y)\ \ p \ satisfies the inequality 



(s) < C 



' f (t){r)dr + C" f E( sup \detV y Z e (y)\ p sup \V y Z 6 (y)\ p ) dr. 

Jo Jo \6£[0,r] 8e[0,r] J 
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By Cauchy-Schwarz inequality and the estimates ([20]) and (|22|) with 2p instead of p, the 
above expectation is seen to be bounded uniformly in r G [0,i],y G M d . We deduce by 
Gronwall's lemma that 

E ( sup \V y [detV y Z s (y)] \ p ) < K i3 (p) 
\*e[o,t] / 

for some constant KiAp) as required, and conclude the third asserted estimate. 

□ 

Proof of Lemma 14. 21 We again consider M = 2, i = 1 and omit the superscript 1 in the 
process Xj t (x), the inverse flow and its derivative. By Lemma 13.91 we can do the change 
of variables Xot{x) = V m the integral with respect to dx inside the expectation in (j!4p . 
Using the semigroup property of the flow and its inverse (see [12]) together with Fubini's 
theorem (thanks to Lemma l4.3p . we deduce that for a.e. y G M , 

£(y) = E[V(t,y)} (23) 

where $?(t,y) is the random function 

(t, y) := exp { jT(rifar,t(l/)) " C 11 * £ far,t(l/)) - C 12 * £ for,t(l/)))*j ^(vo,t(y))detV y r, ,t(y)- 

Notice that we have used the fact that detV y r]Q t t(y) > 0, which follows from detV y rj r! t(y) ^ 
for all r G [0,i] and r i— >• ^yVr,t(y) Dem g continuous and taking the value 1^ at r = £. The 
bound on sup^o^i ||^||oo readily follows from the previous identity, the assumptions on £q 
and the second estimate in Lemma 14.31 

The function y i— > ^(t,y) is moreover of class C 1 by Lemmas 13.91 and 14.31 Since C 11 and 
C 12 have bounded derivatives we deduce that, a.s. for all y G M. d , 

w(t,y) =ex P j jTnCMi/)) - c 11 * £Wy)) - c 12 * e?(^,t(y))dr| 

x eo(^o lt (y))^V7 ?0lt (y) / VV,t(l/) [Vri - (VC 11 ) * £ - (VC 12 ) * e 2 ] (r?r, t (y))* 

JO 

+ vv*(v)v (m,t(y)) det ^voM + toimAv))? [detv V0:t (y)} . 

(24) 

From Lemma 14.31 we get thanks to Cauchy-Schwarz inequality that 

E f |detV»&,t(i/)| J \Vrjr,t(y)\ dr + |V»»,t(y)| |<fe£V»fc,t(y)| + |V [detVrj ,t(y)} |) < oo 

so we can take derivatives under the expectation sign in (|23|) and deduce the existence of 

V^(y)=E[W(t,y)], 
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and moreover that sup te j 0j T] 1 1 I loo < °o. Similarly, sup te r 0T ] ||V£ 2 ||oo < oo. We can now 
rewrite ([24l as 



V*(t,y) =exp U* nfatW) - c 11 * £(vrAy)) - c 12 * £(VrAy))dr\ 

x ^(Vo,t(y))detVvo,t(y) f VV,t(y) [Vrj - C 11 * V# - C 12 * V£ 2 ] (^(y))cZr 

J o 

+ V*r ?0i4 (y)V (r»,t(i/))<fe*Vi» lt (y) + $(fjo,t(!/))V [de*V»fc, t (y)] • 

(25) 

Since ||C^||i = C 1 -?, we have ||C^ * VC2||e» < c^'llV^'Hoo for all s G [0,T]. Taking expecta- 
tion in (|25p and using the estimates in Lemma 14.31 we deduce that for all t 6 [0,T], 

IIV^IU < C" rCHV^Hoo + ||V^||oo)^r + C" 
Jo 

for constants C"',C" > not depending on the functions C lj or doing so only through 
their L 1 norms c lj . Summing the later estimate with the analogous one for ||V£^||oo, we 
conclude thanks to Gronwall's lemma. 

□ 

We are now ready for the 

Proof of Theorem I4.lt Again, we write the proof in the case M = 2. Let e > e > 0. 
To lighten notation, we denote simply by £ = (£ , £ 2 ) and £ = (£ 1 ,£ 2 ) solutions of system 
© in [0, T] respectively with C y = c l3 ip £ and C %3 := c^tpg. Proceeding as in the proof 
of Proposition 13.41 we deduce that for all bounded Lipschitz continuous function (p with 

IMUs < 1, 



ml(pL-PoM 2 +c 



pl t -pl t Mx)e s (x)dx 



+ 



P s >(s)(e a 1 (x)-^(x))dx 



+ 
+ 
+ 
+ 
+ 
+ 



c ii *(e s -e s )(x)PiMx)e s (x)dx 
c 11 * e s (x)p}Mx)(e s (x) - e s (x))dx 



+ 

2 



-i 2 



c 11 * e s (xmMx) - p s L Mx))e s (x)dx 



-i 2 



[C LL -C lL ]*£{x)PiMx)£{x)dx 



1 2 



+ 



c u * e s {xmMx) - p s L Mx))e s (x)dx 



c u * z:(x)p s i Mx)(e s (x) - e s (x))dx 



[c i2 -c i2 ]*e s (x)p s i Mx)e s (x)dx 



>ds. 



(26) 
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Thanks to Lemma EL"2l Pg t <pCl is a bounded Lipschitz function with Lipschitz norm bounded 
independently of e, e and s,t £ [0, T]. We thus can rewrite and bound the first term in the 
third line of of (12611 as follows: 



^ C||£s ~~ £s WjOB* 



for some constant C > not depending on e nor on e. The second term in the third line is 
controlled by 

cwesiLuiwiv sup \piMx) - PiMx)\ 2 < c\\e s \\ue s \\ 2 Tv f 11^-^11^+11^-^11^* 

thanks to Lemma [3.51 b). The term in the fourth line is easily controlled by C\\£~ ~Cl\\'ci3*- 
Using the fact that, by Lemma 14.21 £■? has uniformly bounded derivatives independently of 
e > and s G [0, T], we deduce by the assumption on 7 that 



sup IfC^-C 11 ]*^)!* < C\£-£\ 



and a similar upper bound then follows for the term in the fifth line of (|26|) . The last three 
lines can be bounded in a similar way, and the first line on the right hand side is bounded 
in terms of dual Lipschitz distances by arguments as in Proposition 13.41 Proceeding in a 
similar way as therein, we now obtain the estimate 

i=l,2 

where the constants do not depend on t £ [0, T] nor on e or e. The latter estimate can thus 
be strengthened to 

sup ii£-£ii£b.< cTe su p ne-eii£ B *^ + c|e-ef , 

re[0,t] ■/ i=12 rG[0,s] 

which when summed with the corresponding estimate for i = 2 yields 

E sup lie-CII^<C|e-e-| 2 (27) 

i=l,2^[0,T] 

by Gronwall's lemma. Therefore, as £ goes to 0, the sequence (£ £ ) £ >o is Cauchy in the Polish 
space C([0, T], M. ) and thus converges to some weakly continuous measure flow. Dunford- 
Petits criterion for weak compactness in L , together with the uniform bounds both in L 1 
and L°° for (£f%>o imply that the components of the previous limit have densities for each 
t 6 [0, T] satisfying the same L 1 and L°° bounds, and which we denote by u\(x). 
Weak L 1 -convergence is however not enough to identify u as a solution of (|17p and some 
regularity of the limit will be needed to do so. We denote by Pl t {x,dy) the semigroup 
associated with the measure flow (u l t (x)dx) tG [ ^] as previously and set for cp with ||</?||,C8 < 1, 







^(t,ip) := (ul,ip) - <&P„» - / / (ri(x) - c u ul(x) - c 12 u 2 s {x)) P s >(x)^(^)^- 
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Since £ = £ e satisfies 



(n(x) - C u * - C 12 * £ s 2 (x))P s 1 i ^(x)£ s 1 (x)dxd S = 0, 



proceeding in a similar way as done to obtain the estimate (I26p . we deduce that 



i* 1 ^)! 2 < E su p ne-<n^+ T f / [c li 

i=12 r6[0,t] JO U 



* <{x) - c li <(x)]P s 1 ,99(x)n](x)dx 



ds 



(28) 



(the last term corresponding to the sum of the fifth and last lines in (126p when e = 0). Now, 
by Lemma I3~2| there exists a constant if > (independent of e > ad s G [0, T]) such that 
for any e > one has ) l , d Xl ip)\ < K\\tp\\x for any C°° compactly supported function 
ip. By letting s — > 0, the same bound is satisfied by and by standard results on Sobolev 
spaces (see e.g. Brezis [3] Proposition IX. 3) we get that u has distributional derivatives 
in L°° and that u is Lipschitz continuous with Lipschitz constant less than or equal to K. 
Since C x% * u l s {x) — c ll u\{x) = c u f j(z)[ul(x + ez) — u l s {x)]dz and f j(z)\z\dz < oo, we 
deduce that \\C lt * u\ — c ll u\ ||oo < Ce, which combined with the bound 



V sup Ui-4\\ 2 CB *<CE 2 

following from fl27|) yields |*(i,<^)| < Ce, for all e > 0. That is, u solves f[T7|) . 

Let us finally prove the uniqueness of the solution u. Recall that Pl t (x, dy) denotes the as- 
sociated diffusion semigroup. Consider a second solution v in C([0,T],J\4 M ) with bounded 
Lipschitz norm and associated semigroups denoted P l st {x,dy). Then, 



+ 



+ 



+ 



+ 



{u\ - v},ip) 2 
<(ul(Pl t -Pl t )ip) 2 + C 

c u (ul-v l s )(x)PlMx)ul(x)dx 
c ll v 1 s (x)Pl t ip(x)(u 1 s (x) - v\{x))dx 
c l2 (u 2 s - v 2 s ){x)Pl t ip{x)u]{x)dx 
-} 2 v 2 s {x)Pl t Lp{x)(u](x) - v]{x))dx 



(Pit ~ Pl t Mx)ul(x)dx 



+ 



p lMx)(ul(x) - v\(x))dx 



+ 



c n vl(x)(PlMx) - P}Mx))ul{x)dx 



+ 



-12 2 



v 2 s (x)(PlMx) - Pl t <p(x))ul(x)dx 



>ds. 



(29) 



The functions P} t ipu\ and Pl t ipv\ having Lipschitz norm bounded independently of s, t € 
[0,T], the first term in the third line and the term in the forth line above are bounded by 
C\\u~ — fgHjcg*. The second term in the third line is bounded by 



ClksllLhslli SU P 



x-eit 



PlMx) - PlMx) <c\\ 



.1 M 2 ||. 1 1 1 2 
ll 



.1 „.1||2 



u^ — v. 



r\\CB' 



+ \\U r — V. 



2 1 1 2 
r\\CB 



*dr 
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by Lemma l3.5l b). The last threes lines are similarly dealt with, and we can easily conclude 
as in Proposition 13. 4[ 

□ 
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